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ABSTRACT 

We analyse the quark mass hierarchy and CKM matrix using the universal Yukawa 
coupling model with small violations precisely. We estimate the ranges of the values of 8 
violation parameters (5", 62, S^, 6f, 62, S^, 02, 03) in our quark mass matrices satisfying 
quark mass ratios and CKM matrix, where 02, 03 are phases. Without these phases, the 
solution satisfying quark mass ratios and CKM matrix is not obtained. These parameters 
obtained can explain the CP violation effects. 
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1. Introduction 



The origin of the mass hierarchy of quarks and leptons has been investigated by the 
various theories beyond the standard model (SM) by many authors Although it 

is necessary to study precisely the theories beyond the SM, in the present circumstances 
where there are the precise analysis of — mixing , the CP- violating parameter e of 
the — system and the determination of the top-quark mass, one should analyse the 
mass hierarchy of quarks and Cabibbo-Kobayashi-Maskawa (CKM) matrix independently 
of model assumed. Furthermore, many attempts beyond SM |^ constructing mass matrix 
patterns at the GUT scale in SUSY theories or at the string scale in string models, 
although quite successful, cannot produce results in complete agreement with precise low 
energy data. Thus, for the model building beyond SM, the analysis using only the minimal 
qualities to be sure at present is now to be very necessary. 

For quark mass matrix patterns in low energy, there are Fritzsch type, Stech type 
model the democratic type model and the universal Yukawa coupling type model 
0]. We adopt a quark mass matrix like the democratic and the universal Yukawa 
coupling type models 0] with the small violations from the universality which cause the 
mass hierarchy. Our model does not take any assumptions on the violations and treats 
violation parameters as free parameters. First, we see the mass hierarchy mechanism 
in the limit of the universal Yukawa coupling. The (m, c, t) and (c?, s, h) quark mass 
matrices are expressed, under the universality of Yukawa coupling strength, as 



where F** and F'^ are real, and are not assumed universal. This F" and F'^ unuiversality is 



quarks acquire their masses through the couplings to two different Higgs multiplets. It 
is well known that this type of the mass matrix is diagonalized as diag[0, 0, 3F^] by the 




(1) 



guaranteed by e.g., a minimal supersymmetric gauge model |0 in which the up and down 
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orthogonal matrix Tq; diag[0, 0, SF^] = TqMTq ^, where Tq is 

(2) 



Tn 



1 1 2_ 

\ ^ A. ^ , 

\ V3 / 



Thus the type of Eq. (1) gives the large mass gap between the heaviest quark and other 
two quarks. 

Next, we introduce small violations with phases of Yukawa coupling strength in uni- 
versal coupling (Eq. (1)) as 



/ 1 l-6te'^i l-5|e'^2 

1 - 6le-'-^' 1 1 - S^e'^s I , {q = u, d) (3) 

V 1 - 6le-'^'^ 1 - 5fe-^^3 1 



where 5| are small real violation parameters 

5t'' « 1. (z = 1, 2, 3) (4) 

We do not take any assumptions on the violation parameters except that (5|"s are very 

small. Introducing small violations of coupling, we assumed that the violations are caused 

from the coupling between different quarks. Here it should be stressed that the large mass 

differences are produced by the universal coupling then what distinguishes the quarks is 

not the masses but other characters. Thus the assumption that the diagonal elements 

of couplings between same quarks are same each other and do not have violations is 

reasonable. In the parametrization (Eq. (3)), we used the — sign before 61 because, in 

this notation, 6f are allowed only to be positive in the quark mass analysis as shown later. 

Branco, Silva-Marcos and Rebelo 0] studied the type of mass matrix (3) but they equate 

the type of this mass matrix to the quark mass squared. 

Under the assumption of small violation in universal coupling, we can get the second 

mass gap between two degenerate zero mass states which are taken from the universal 

mass matrix (1). This is shown from the mass matrix (3) neglecting phases 

1 1-61 1-6^2 
M^ = T''\ 1-61 1 1 - -^s I • iQ = u, d) (5) 
1-6', 1-61 1 



This mass matrix is transformed by Tq to 
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6t iSl-6l)/VS {-5I + 5I)/V6 

(5^ - Sl)/V3 {-61 + 251 + 25l)/3 {-251 + + ^l)/^V2 
{-51 + 51)/VQ {-25l + 5l + 5l)/?,V2 {<^ - 251 - 251 - 251) / ?, 



) 



(6) 



then if 5? < 5| ^ 5^ < 1, three eigenvalues turn to (0, 45|rV3, (3 - 45|/3)r«), approxi- 
mately. This tendency is certified and the allowed ranges of {51, 52, 5^) are determined 
precisely from the analysis of quark mass ratios and CKM matrix in the following numer- 
ical study. 

In last section (3. Discussions), we will comment on the difference between the results 
of our model and others 

2. Numerical analysis 

The mass matrices (3) contain 6 violation parameters for {u, c, t) and {d, s, h) sector 
except for the V^, respectively. First, we consider 3 parameters {51, 5^, 51) case (Eq. (5)) 
for simplicity and later we will consider the case containing phases. We diagonalize the 
mass matrices (5) to the diag[m^i, m^, m^] and diagfm^, m^, m;,] for q = u and d by the 
unitary matrices T(5", 52, 5^) and T{5f, 52, 5^), respectively. 



Eigenvalues of the mass matrices are not the physical masses but the parameters in 
the Lagrangian. Theses quark masses (eigenvalues) are running masses which should be 
all taken on a single energy scale. In order to estimate the parameters {51, 51? '^3)5 we 
use the quark mass ratios. These mass ratios are, to a good approximation, independent 
of the energy scale, then the scale can be arbitrarily chosen §]. 

For the values of the light and medium heavy quark u, d, s and c masses, we use the 



T{51, 51 5l)M^T-\5l 51 51) = Ml 



{q = u, d) 



= diag[m„, m^, m^], = diagfrn^, m^, nib]. 



(7) 
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world average cited in Ref. 0]; 

m„ = 2 - 8 MeV, = 1.0 - 1.6 GeV, 

mrf = 5 - 15 MeV, m, = 100 - 300 MeV, (8) 

because these values are almost similar to the mass values at the scale fi = IGeV P|; 
ruu = 5.1 ± 1.5MeV, = 8.9 ± 2.6MeV, = 175 ± 55MeV, = 1.35 ± 0.05GeV. 

For the heavy quark b and t masses, we estimate the running mass mq{ii = IGeV) 
related to the physical mass ■mP^y'^ in the first order QCD as 

mfy' = mq{fx = mg) 



4 

1 + —asifi = ruq) 
on 



(9) 



where the running coupling constant as{fi) and the running mass rrig^fi) are expressed as 

gifif 4n ( (3i \nL\ 



Air (3oL \ (3i L 



_( 2/?i7olnL + l 871 n 

/3o = ll-^iV/, 70 = 2, = 102 -^Nf, 

102 5 ^ u2 

7i = jVf, L = in—. 

12 18 ^' A2 

The estimated rriblfi = IGeV) and rritlfi = IGeV) from the physical mass m 

4.3 ± 0.2 GeV and mf^' = 174 ± g GeV § are 



phys 
b 



A = 0. IGeV A = 0.2GeV 
mbifi = IGeV) = 5.08 ± 0.28 GeV, 5.49 ± 0.29 GeV, 
mtifx = IGeV) = 289 ± 41 GeV, 327 ± 48 GeV, (11) 

for the flavor number Nj = 3. We write the values of mb(/i = IGeV) and mt{fi = IGeV) 
for the renormalization group invariant scale A = O.lGeV and A = 0.2GeV cases, because 
these mass values are sensitive to the values of A. Hereafter, we write the mass mq{fi = 
IGeV) as rriq. From these mass values, we get the quark mass ratios, 

— = 0.0038 ±0.0025, — = 0.050 ± 0.035, 

rrir m. 



— = 0.0042 ±0.0013, 



— = 0.038 ±0.019, 



(12) 



where we used the average values of m5(/i = IGeV) and mt(/i = IGeV) for A = O.lGeV 
and A — 0.2GeV and involved the deviation from the average value in errors. 

We estimated numerically the allowed regions of {Sf, 51, S^) satisfying the constraint 
in which the ratios of the eigenvalues (m„, ttIc, mt) and (m^, nis, mi) of the mass matrices 
(5) are included in the experimental ranges of quark mass ratios (12). We showed the 
allowed regions for (m„, m^ mt) sector in Figs. 1(a), (b) ,(c) and for (m^, m^, m;,) sector 



The Figs. 1(a), (b), (c) represent the allowed regions of {S'!^-,^'^) plane corresponding 
to the S'l - 0.00005, 0.0001, 0.0004 for {u, c, i)sector, respectively and Figs. 1(d), (e), 
(f) the allowed regions of (^2, 5f) plane corresponding to the 5f — 0.005, 0.01, 0.02 for 
(ci, s, h) sector, respectively. For 5" < 0.000012 and 5f < 0.00085, the allowed regions 
for (^2 ) ^3) ^iid (<^2) ^i) plane do not exist, respectively. It is seen in this Fig. 1 that the 
allowed regions for (51, ^f) are symmetric with respect to the interchange between 51 and 
6^ . This symmetry is found easily in the approximate expressions for the eigenvalues and 
eigenvectors of the mass matrix (5). The eigenvalues are 



in Figs. 1(d), (e), (f). 



Fig. 1 (a), (b), (c), (d), (e), (f) 




list + 61 + 61)-^^"] 

list + si + 5i) + l^^'\ 

3-^(5?±5|±5I)] 



(13) 



where 



(25? - 5| - Sir + 3(5| - 5|) 



2lV2 



(14) 
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and the corresponding eigenvectors are Uf, U2, 

THSl5l5l) = [iU!),iUl),iUl)], 

cos 01 sine" A*? cos ^« + /i« sin ^« 

T{Sf,6l5l)^\ -sin^« cos^-^ -A^sin^^ + Zi^cos^-^ | Tq, (15) 

-X" 1 



where 



1 



A" = 77^(^1 - SI), /i'' = —=(25? - 51 - 51), 



^^-Itan-^S^i^ (16) 



Though these expressions are obtained approximately, the allowed regions for (^2, ^3) 
obtained from these approximate expressions are almost same as those in Fig. 1. 
Next we consider the CKM matrix V, 

V = T{5l,5^2,5l)T\5i,5i,5i). (17) 

The matrix elements of V are determined by various experiments, for example, nuclear 
beta decays, K^.^ decays, neutrino and antineutrino production of charm off valence d 
quarks, neutrino production of charm, semileptonic decays of B mesons produced on the 
T(4S') hh resonance and etc. The absolute values for these matrix elements are tabulated 
as 

/ 0.9747- 0.9759 0.218 - 0.224 0.002 - 0.005 \ 
yexp ^ Q _ 0.224 0.9738 - 0.9752 0.032 - 0.048 . (18) 
V 0.004- 0.015 0.030 - 0.048 0.9988 - 0.9995 / 

We calculated numerically the allowed regions of {51, 52, 5^) satisfying the restriction in 
which the absolute values of matrix elements of V in Eq. (17) are included in the exper- 
imental range of matrix elements (Eq. (18)). First, we estimated the allowed regions of 
{5i,52,5l) independent of the experimental constraint of ranges of the mass ratio (12). 
We showed the allowed region of (^2, ^3) for {u, c, t) sector in Fig. 2(a) and of (5|, 5|) 
for {d, s, b) sector in Fig. 2(b) fixing the 5" and 5f as 5^ = 0.0001 and 5f = 0.01. For 
other values of 5"^ and 5f, 5"^ = 0.00005, 0.0004 and 5f = 0.005, 0.02, the allowed re- 
gions for {52, 5^) and (^j, 5^) are almost similar to those of the case 5" = 0.0001 and 



5f — 0.01 shown in Fig .2. All combinations of all points in allowed region of {62, S^) with 
all points in allowed regions of , 5^) are not allowed but the restricted combinations 
between some points in allowed region of {82, 5^) and some points in allowed region of 
(^2) ^3) are allowed. For example, the combinations of a point in (^2, 5") plain (shown 
by a large dot) with points in the area shown by the large dotts in {62, 5|) plane are 
allowed. In Fig. 2, though we showed the solution corresponding to the case, 62 > (^3 and 
S2 < 5f , there exist also the solutions corresponding to the cases, 6^ < ^3 and 62 > S^- 



Fig. 2 (a), (b) 



As shown in Fig. 1 and Fig. 2, there is no common regions satisfying both constraints 
of the mass ratios and the CKM matrix in (^2 , 5f ) plane of {d, s, b) sector. This fact is 
easily understood from the analytic expressions for V. Using the approximate expression 
(15) ior T {61,82, 51), we can get the approximate expression for V 

/ cos(r - 6'^) sin(e" - 6'^) 

V ^ - sin(^" - 9"^) cos(^" - 0^^) 

^ -(A" - A*^) cos 9"^ - (//" - //<^) sin 9'^ (A" - A'') sin 9"^ - (//" - //<^) cos 9^^ 

(A" - A"^) cos + (/i" - /) sin ^ \ 

-(A" - A'^) sin + (/x" - /) cos 1 . (19) 

From this expression, for ^2 ~ (^f » (^1, we can get the ratio iKb/V^bl ~ | — (A" — 
A'^)sin6'" + (/i" — ii'^')cos6"\ ~ and then 62 ^ > 0.20 from the experimental 
ratio \V^^/V^l^^\ > 0.032. On the other hand, from the ratio of mass eigenvalues (13), 
for 5f fa 5f > 5f, we can get m^/mf ^ 2(5f + (5f)/9 then 5f 5f < 0.13 from the 
experimental range ms/rrii) = 0.038 ± 0.019. We comment here on the difference between 
the values calculated by Eq. (19) and by numerically exact procedure. The values for 
Vus and Vcd elements calculated by Eq. (19) are different from the values calculated by 



numerical and exact procedure about 2%, and for Vub, Vcb, Vta and Vts about 20% and for 
Vud, Vcs and Vtb about 0.1%. 

Because of the fact that there is no common region in [82, 5f) plane satisfying the 
quark mass ratios and the CKM matrix, we consider the case containing the phases ifl in 
the quark mass matrix as Eq. (3). Although there are 6 degrees of freedom for phases 
only two phases and ip'^ are considered in our analysis because only 5f and 5| in the 
violation parameters are about 0.1 and other parameters are extremely small {5f ~ 0.01, 
5" ~ 0.0001, ^2,3 ~ 0.01), then the phases with these other parameters scarcely contribute 
to CKM matrix in contrast to two phases and . We parametrize the (d, s, h) sector 
quark mass matrix using the very small phases 02 and 03 instead of the phases (^2 ^^'^ 



I 1 

\ (1 - 



l-5f 
1 



:i 



-i<f>3 



(1 - 5^)e''f'^ \ 
(1 - (5f)e^'^3 

1 



<Pi « 1, {i = 2, 3) (20) 



J 



and for M" we use the type of Eq. (5) with no phase. The approximate expressions for 

the eigenvalues of the mass matrix (20) are the same expressions for mf, m2 and m'l 
as Eq. (13), but the expression for is changed to containing the phases 02 and 03 as 
follows. 



1/2 



e = [{26t - - Sty + 3(^2' - ^3)' + 3(02 - 03)' 

The expression for CKM matrix is given in this approximation as 



V = T{5lSlS^)T\St,SiSi<P2,<p3) 

j cos ^"c'^ + sin ^^s'^* - cos ^"s"^ + sin O'^d^ 

fti -sin^V + cos^^s'^* sin ^"5^^ + cos ^"c'^ 

V -(A" - A'^*)c'^ - (;u" - ^'^*)s'^* (A" - A'^*)s'^ - (/x" - /x'^*)c'^ 

(A" - A"^) cos ^" + (/i" - /i'^) sin ^« 
-(A" - A'^) sin + (/i" - /) cos ^" 
1 



(21) 



(22) 



where 



■V^{(52-'^l)-^(02-03)} 



/ = ^ {2St - {S', + St) - 3z{h + 03)} , (23) 

and A" and are given in Eq. (16). We calculated the allowed regions of {5^, 82, (^3) 
and {5f, 62, 62, 02, ^3) numerically and exactly satisfying two constraints of the ex- 
perimental values of ranges of the quark mass ratios Eq. (12) and the CKM matrix V^^^ 
Eq. (18) and showed these regions in Fig. 3. We showed the case ((^", 5f, 02, 03)=(0. 00005, 
0.005, -4°, -4°) in Fig. 3(a), (0.0001, 0.01, -4°, -4°) in Fig. 3(b). Solutions correspond- 
ing to the cases other than (02, 03) ~ (—4°, —3° ~ —4°) do not exist. 



Fig. 3 (a), (b) 



In order to see the effects of the CP violation, we rephase the CKM matrix V — 
T{5'1, 5^, 5f)T\5f, 5^, 5f , 02, 03) to the standard parametrized CKM matrix where the 
matrix elements V!^, and Vj^ are real number, by using the rephasing matrix 

and Pd as 

T/« = PuVPI 

Pu = diag[e*"', 1, e'^'], Pa = diag^, 1, (24) 
The parameters a, a', (3 and /3' are determined as 

a = tan „ ^, tan — — — , a = — tan 



^ _i ImV^fe , _i ImV^ft _i ImVtf, 

/3 = tan — — , /3 =tan ——-tan — — . 



(25) 
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In the standard parametrized CKM matrix in which the element V^^ is almost real as 
recognized in the Wolfenstein parametrization 0, the parameters p and t] characterizing 
the CP violation which are the vertex coordinate of unitarity triangle are expressed as 

^_MYrv^ „ Im(KfKf.) ,,,, 

The phenomenological constraints for parameters p and rj has been examined by Pich and 
Prades ^ using the recent information on the non-perturbative hadronic inputs needed 
in the analysis of — mixing and the CP-violating parameter e of — system. 
They gave the results of parameters p and r] for the best estimate set of input parameters. 
We showed our resuhs of (p, r]) for {5^, 6f, 02, h) fixed as (0.00005, 0.005, -4°, -4°) 
and (0.0001, 0.01, —4°, —4°) in Fig. 4, besides the Pich and Prades results which 
are surround by circles centered at (0, 0) and (1, 0) and hyperbola correspond to the 
input parameters {xd = 0.76 ± 0.06, mf'" = 174 ± 16 GeV, r(5^) = 1.61 ± 0.08 ps, 
T{B^a)W = (3.9 ±0.6) X 10^GeV-\ \Vub/Vcb\ = 0.08 ±0.03, Bk = 0.50 ±0.15, Cb/A = 
2.0 ± 0.5) [Q. From this Fig. 4, we can say that the values of parameters (5", 62, S^) 
and {6f, 02, ^3) in allowed regions shown in Fig. 3 are almost consistent with the 

phenomenological CP violation results. 



Fie. 4 



We summarized the ranges of parameters (5", 5" = (^2 ± ^3)/'^, 6^ = 62 — ^3) and 
{5f, 5^ = {5^ + 5i)/2, 5i = 5!j-5l 0+ = (02±03)/2, 0_ = 02-03) obtained in previous 
analysis in Table 1. The compound signs with the values of 5", and 0_ correspond to 
each other. Of course, all combinations of values in the ranges shown in Table 1 are not 
the solutions but special combinations are the solutions. 
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Table 1: The ranges of parameters (5", 5", 5", Sf, S"^, 51, 0_) satisfying the mass 
ratios (Eq. (12)) and CKM matrix (Eq. (18)). 



uct sector 


dsb sector 


61 = d^- 


0.00001 ~ 0.0004 
0.0064 ~ 0.0125 
±(0.0 - 0.0043) 


si = 5i - 5i 

4)% = (4)2 + 4)?)I2 

<?'L = 02 - Q:i 


0.001 ~ 0.015 
0.040 ~ 0.129 

±(-0.038 0.006) 

-4° ~ -3° 

±(-r - 0°) 



We show the typical solutions, and the mass ratios, CKM matrix elements and (p, rj) 
corresponding to these solutions; 



solution A 



solution B 



6^ = 0.00005, 6^ = 0.01, 6^ = 0.009, 

5f = 0.005, 5^ = 0.054, 5f = 0.08, (/)2 = -4°, 03 = -4°, 



m,, rrir 
— = 0.0019, — 
rrir 



0.0042, ^ 



0.025, ^ 
nib 



0.031, 



mt 

f 0.9755 0.2198 0.0036 
V= 0.2196 0.9750 0.0340 

V 0.0077 0.0333 0.9994 
p = 0.088, T] = 0.47, 

j S^ = 0.0001, = 0.01, = 0.009, 

\ df^ 0.01, 5^ = 0.07, Si = 0.102, 02 = -4°, 03 = -4°, 



^ = 0.0058, — 
rric rut 



Try J Tfi 

0.0042, — = 0.058, — = 0.040, 
ms rUb 



/ 0.9753 0.2210 0.0043 

V= 0.2210 0.9747 0.0347 

V 0.0087 0.0339 0.9994 

p = 0.022, r? = 0.56. 



(27) 



3. Discussions 

We analysed precisely the mass hierarchy of quarks and CKM matrix in the law 
energy by the universal Yukawa coupling model with small violations (Eq.(3)). Violation 
parameters estimated are 8: {Si, S2, S^, Sf, S2, Si, 02, 03), and estimated values of these 
are tabulated in Table 1. Other phases than 02 and 03 do not contribute to our present 
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analysis because of extreme smallness of (5", 62, S^, 6f). We fitted 8 violation parameters 
to 8 experimental data: 4 quark mass ratios (Eq. (12)), 3 mixing angles determined by 
CKM matrix elements (Eq. (18)) and 1 phase determined by CP violation which relates 
to p and 1]. 

Here we comment on the differences between our model and others JH-Q]. The mass 
matrices depending on models adopted and CKM matrix are connected through the uni- 
tary matrices Tu^d depending on models as follows; 

TuM'^T-' = M^ = diagK, m„ m^], 
TdM'^T^^ = Mi = diag[mrf, m„ m^], 

V = T^Tl (28) 

From V = TuT\ and V^^^ ~ 1, Fritzsch type , Stech type model and many other models |ll|], 
i adopt the unitary matrices Tlf^'"''' as Tlf^'"'' ^ 1 and then (M")^"*^^^'^ ^^p"^ = 
(j^Yitzsch type)-i^jag[^^^ ^^jj^Fritzsch type ^ diag[m„, m^, m^] and (M'^')^"*^"^^ *yp*^ = 
^jiFntzsch type-j_i^j^g|-^^^ mbjTj"*^'''^'^ ^■^P'' ^ diag[mrf, m^, mi]. However, if we change 
Tu — ^ TuS and Td TdS where S is some arbitrary unitary matrix, the CKM matrix V 
remains unchanged. The democratic model and the universal Yukawa coupling model |0| , 
1^, in fact, use the following uitary matrices T^^^]^"^^"^^ coupling type 

Tiuniversal coupling type rj-iFritzsch typerj-i lnc\\ 

where Tq is the unitary matrix defined in Eq. (2). CKM matrix does not depend on the 
unitary matrices T^^d adopted but the weak interaction eigenstates do on it then on the 
mass matrices adopted. 

Even in Fritzsch type models, there are many parametrizations and all parametriza- 
tions cannot explain precisely the present quark mass hierarchy and CKM matrix. For 
example, original Fritzsch model can not explain the observed large top quark mass and 
modified model (see the literature of B.Dutta and S. Nandi in [M) can explain the large 
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top quark mass but has to use the up quark mass matrix in which the (2, 3) and the 
(3,2) elements are unequal Recent Peccei and Wang analysis (the hterature in [^) uses 
the mass matrix M"^ = T~^{diag mf[^utA^, ^ct^^, ~ diag mt[^ut^'^, xictX^, 1] and 

= T7^{diag m,[aA^ ^.^A^, l]}Trf ^ diag m,[aA^ 1], where U = 0-49, ^ct = 

1.46, ^db = 0.58, ^sb = 0.55 and A is the parameter in Wolfenstein parametrization of 
CKM matrix, and takes the mass matrices explaining the law energy data precisely. But 
the principle to take their mass matrices is not so clear. 

We will analyse the problem of neutrino mixing using the present our model in next 
work. As we mentioned above, the weak interaction eigenstates depend on the unitary 
matrices Tu^d adopted, then the analysis of mass matrix involving the lepton sector like 
the neutrino mixing problem will give the clue to check the validity of our model and 
other models. 
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Figuare captions 

Fig. 1. The allowed regions for (5", ^3) i^i^ ^2) '^3) satisfying the mass ratios 
(Eq. (12)). (a), (b), (c): The allowed regions of (525*^3) plane corresponding to the 
5" = 0.00005, 0.0001, 0.0004 for {u,c,t) sectors, respectively, (d), (e), (f): The allowed 
regions of (^2, 5|) plane corresponding to the 6f = 0.005, 0.01, 0.02 for (d, s, h) sectors, 
respectively. 

Fig. 2. The allowed regions of (^J*, 82, ^3) and {5f, 62, S^) satisfying the experimental 
CKM matrix elements (Eq. (18)). (a): The allowed region of {62, ^3) plane corresponding 
to the 6"^ = 0.0001, 6f = 0.01. (b): The allowed regions of {62, S^) plane corresponding to 
the 5i = 0.0001, 6f = 0.01. All combinations of all points in allowed region of {62, S^) with 
all points in allowed regions of {62, S^) are not allowed. For example, the combinations 
of a point in (^2, S3) plain (shown by a large dot) with points in the area shown by the 
large dotts in (5|, 5|) plane are allowed. 

Fig. 3. The allowed regions of (5", 62, S^) and {Sf, S^) satisfying the mass ratios 
(Eq. (12)) and the experimental CKM matrix elements (Eq. (18)). (a): the allowed regions 
of {6^, 6^) and {6^, 6^) planes corresponding to 6"^ = 0.00005, 6f = 0.005, (j)2 = -4°, 03 = 
—4°. (b): the allowed regions of (^2 , ^3) and (^2 , ^3) planes corresponding to 5" = 0.0001, 
6f = 0.01, 02 = -4°, 03 = -4°. 

Fig. 4. The parameters (p, rj) satisfying the mass ratios (Eq. (12)) and the experimental 
CKM matrix elements (Eq. (18)) for the (5^ 5f, 02, 03) fixed as (0.00005, 0.005, -4°, -4°) 
and (0.0001, 0.01, —4°, —4°). Area surrounded by the circles centered at (0, 0) and (1, 0) 
and hyperbola is the allowed region for CP violation given by Pich and Prades [§. 
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